Nonlinear effects and the behavior of total hadronic and photonic cross sections 
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In this work we use an eikonalized minijet model where the effects of the first nonlinear corrections 
to the DGLAP equations are taken into account. The contributions coming from gluon recombi- 
nation effects are included in the DGLAP+GLRMQ approach for the free proton in the context 
of saturation models. The parameters of the model are fixed to fit total pp and pp cross sections, 
including the very recent data from LHC, HiRes and P. Auger Collaboration. Glauber and multiple 
scattering approximations are then used to describe the inclusive inelastic proton-air cross section. 
Photoproduction cross sections, without change of parameters fixed before, are also obtained from 
the model using Vector Meson Dominance and the Additive Quark Model. We show and discuss the 
main results of this study as well as the implications of saturation effects in the behavior of total 
hadronic and photonic cross sections. 

PACS numbers: 11.80.Fv, 24.85. +p, 25.75.Bh, 13.85.Tp 



I. INTRODUCTION 

The growth of total hadronic cross sections with energy 
has been studied for decades and many phenomenological 
and theoretical efforts have been made to explain it. 

One of the most important explanations for this behav- 
ior was proposed in the 70's |l| based on QCD (Quan- 
tum Chromodynamics): the behavior of hadronic cross 
sections with the center of mass energy (y/s) was very 
similar to the production of jets, indicating that partons 
would be playing a role in these interactions. This basic 
idea has led to an approach, called minijet model, which 
takes into account that the total hadronic cross sections 
can be decomposed as oq + u p qcd , where character- 
izes a nonperturbative contribution to the process (gen- 
erally taken as energy-independent at high energies) and 
&pQCD represents the semihard contributions, calculated 
in perturbative QCD (pQCD) with use of an arbitrary 
cutoff at low transverse momenta pr min (> ^-qcd) [2j|- 

This simple model, however, violates the unitarity of 
the iS-matrix for these process and, consequently, the 
Froissart bound, which states that total hadronic cross 



sections cannot grow faster than In (s) as s — > oo 
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At high energies, the perturbative component of this 
model is dominated by gluons with very small (Bjorken- 
x) fractional momentum and pQCD calculations, based 
on linear QCD evolution (DGLAP Q and BFKL equa- 
tions @), show that the minijet cross-section grows very 
rapidly, dictated by a power-like energy behavior. 

To explain the experimental data quantitatively, this 
idea was reformulated based on the eikonal representa- 
tion to ensure unitarity and to contain this strong level of 
growth. Since then several models (Eikonalized Minijet 
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Models - EMM) have been proposed [6j and, in general, 
good descriptions of the experimental data have been 
obtained. Nevertheless, many questions about the dy- 
namics of interaction between partons in this regime still 
remain open. 

Parallel to these developments, significant progress 
were achieved in theoretical physics of small- a; and the 
results from HERA, with kinematical ranges extended 
upwards in Q 2 (the four-momentum transfer to the pro- 
ton) and also downwards in x, changed our view of 
the structure of the proton Deep inelastic scatter- 
ing experiments showed a rapid increase in the density 
of gluons as x decreases and reinforced the hypothesis 
that this growth could be related to nonlinear effects in 
gluon evolution equations. One now know that for values 
x < 10~ 4 , gluons dominate the hadron wave functions 
but, of course, it is expected that the growth of gluon 
densities "saturates" at a given time. 

The understanding of this expectation is related to the 
momentum transfer, fej_, and, therefore, to the transverse 
size of a gluon (oc l/fc±) in semihard interactions. For 
large momentum transfer, the BFKL evolution predicts 
a large number of small size gluons per unit of rapidity 
produced through g — > gg interactions. For small mo- 
mentum transfer, on the other hand, the produced glu- 
ons overlap themselves in the transverse area and fusion 
processes, gg — > g, become also important. 

This simple scenario shows that a typical scale, Q s , 
called "saturation scale", tell us that these latter pro- 
cesses are small for k 2 j_ > Q 2 S . For high enough momen- 
tum transfer, k\ < Q 2 , however, Q s tell us that the 
recombination of gluons (fusion processes) cannot be ne- 
glected because the gluon density is large and grows with 
lowering x. At very high energies, smaller and smaller 
values of Bjorken-x can be accessed and, under these con- 
ditions, the recombination mechanism becomes more and 
more effective resulting in a decrease in the population 
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of gluons and, therefore, in the idea of "saturation" of 
partonic distributions mentioned above. 

Many studies were made in the latest two decades ex- 
ploring this subject and, currently, one believes that an 
effective theory, the Color Glass Condensate (CGC) 0- 
rToj ] , correctly describes the behavior of very small- a: glu- 
ons in hadronic wave functions by an infinite hierarchy of 
coupled evolution equations for the correlators of Wilson 
lines. 

According to this picture, the (highly dense) system 
formed at these extreme conditions is characterized by 
the limitation on the maximum phase-space parton den- 
sity that can be reached in the hadron wave function 
and an x— or energy-dependent momentum scale, Q s (x), 
which separates dense and dilute regimes. In the low den- 
sity regime the formalism reproduces the BFKL dynam- 
ics for partons with transverse momentum much larger 
than this "saturation momentum" . On the opposite side, 
the saturation scale becomes large for Q s {x) 3> Aqcd 
and the formalism predicts that partons saturate in the 
hadron wave function with occupation numbers of order 
l/a s (Q s ). In this case the coupling constant becomes 
weak (a s (Q s ) 1) and the high energy limit of QCD 
can be studied using weak coupling techniques [TT|, [l2| ■ 

In this work we study the influence of gluon recombina- 
tion process assuming that such a system may be formed 
in hadronic and photonic collisions at very high energies. 
As mentioned above, at very small x-values, parton dis- 
tribution functions are governed by BFKL dynamics and 
this mechanism leads to nonlinear power corrections to 
the DGLAP evolution equations. We also adopt here the 
first nonlinear (GLRMQ) terms calculated by Gribov, 
Levin and Ryskin |13j . and after by Mueller and Qiu 
14], to describe experimental cross sections and make 
predictions. 

In what follows we briefly present the standard formu- 
lations of the eikonalized minijet model for hadronic and 
photonic cross sections, of the inelastic proton-air cross 
section and describe the main ingredients used in our cal- 
culations. We then present the strategy used to fix the 
parameters of the model, show and discuss our main re- 
sults and, in the last section, we outline our conclusions. 



II. EIKONALIZED MINIJET MODEL WITH 
SATURATION EFFECTS 

One of the most important contributions to predict 
the behavior of hadronic cross sections with the energy 
from the QCD parton model was proposed by Durand 
and Pi in the late 80 's using a formalism consistent 
with unitarity constraints. Many QCD-inspired models 
used today have their origins based on this "eikonal" for- 
mulation, which provides a framework where the minijet 
cross-sections are unitarized via multiple scattering. 

In this work we have used an unitarized version of the 
minijet model j Hi where the total, elastic and inelastic 



pp(p) cross sections are given by: 

aZ iP \s) = 2 / d 2 b{\ - e- Im ^cos[Re X (b, s)]} , 



^f (p_) (s)= / d 2 b\l~e l ^\ 2 



«EP« = d%[l-e 



-2 Im x(6,s)] 



(i) 



The eikonal function %(&, s) in the above expressions 
contains the energy and the transverse momentum de- 
pendence of matter distribution in the colliding particles 
and, through the impact parameter distribution in the b- 
space, it is given by s) = Re [x(b, s)] + i Im [ X (b, s)]. 

The real part of X (b, s) represents only about 4% in the 
ratio of the real to the imaginary part of the forward elas- 
tic amplitude for pp(p) processes and therefore, as a first 
approximation, we assume i?ex(6, s) = in this work. 
We also assume that multiple partonic interactions are 
Poisson distributed with an average number separated 
in soft and hard processes in a given inelastic collision, 
n(b, s) = 2Im X (b, s) = n so f t (b, s) + n hard (b, s), and can 
be factorized in b and s as jla |: 



n(M) =W(b, fx soft )a soft (s) 



w{b,^ hard )o h k r d {s), 



k.l 



(2) 



where Wib, /i so /t) and Wib, Hhard)y which represent the 
effective overlap functions of the nucleons at impact 
parameter b, are related to the nucleon form factor 
in hadronic and partonic levels (normalized so that 
JW(b,fx)cPb = 1), and a soft {s) and cr^ rd (s) represent 
the behavior of the total cross sections with energy in 
soft and hard (minijet production (mj)) regimes inpp(p) 
collisions. 

At low energies the hard contribution to the eikonal 
function is very small. In order to describe pp and pp 
scattering at ISR energies we parameterize the soft con- 
tribution as 1171: 



W{b,(i sof t) = -^{fJ, so f t b) 3 K s (}j, soft b), 
»07r 
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where Hsoft is an adjustable parameter, K is the modified 
Bessel functions and is the proton energy in the 
laboratory system (cross-sections are in mb). 

The hard contribution is described by the minijet pro- 
duction in leading order (LO) pQCD where partons are 
produced back-to-back in the transverse plane according 
to the differential cross section UM: 
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da ki 
dy 
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where h\ and I12 denote the colliding hadrons and 
da l ^ kl jdt the subprocess cross sections [19j . 

The rapidities of the final state partons k and I are 
labeled by y (= y\) and 7/2 and the transverse momen- 
tum of each parton by pr (> PT min , the smallest trans- 
verse momentum of parton scatterings allowed). The 
fractional momenta of the colliding partons i and j are 
X1.2 = PT/y/s(e ±y + e ±y2 ), i.e., the incoming partons 
are collinear with the beams. The factor 1/(1 + Ski) is a 
statistical factor for identical particles in the final state. 

In our calculations we have assumed k = 1 and only 
considered the process gg — > gg, gq{q) — > gq{q) and gg — » 
qq (q = u 7 d 7 s). We also have parameterized the "hard" 
overlap functions in impact parameter space, W(b, p gg ), 
W(b,fj, gq = y/^qq^gg) and W(b, p qq ) , as Fourier trans- 
forms of a dipole form factor (see eq. ((3])) [2(J. The (free) 
parameters p qq and p gg represent masses which describe 
the "area" occupied by quarks and gluons, respectively, 
in the colliding protons. 

As discussed before, even in conventional eikonalizcd 
minijet models the rise of the total pp(p) cross section 
with energy is related to the increasing probability of 
perturbative small- a; gluon-gluon collisions: gluon dis- 
tribution functions governed by DGLAP evolution and 
contributions of partons with px > PT min dominate the 
integrand of eq. ((4]) increasing very fast the rise of total 
cross-section with energy. The numerical evaluation of 
this partonic contribution strongly depends upon pr min , 
the chosen set of parton densities (fij/tH 2 ( x i,2,Q 2 )) and, 
of course, their evolution in this regime. 

The main ingredient of our model is the introduction 
of nonlinear terms in the evolution of parton densities 
above. In the context of saturation models, we adopt 
these corrections and make use of EHKQS parton distri- 
bution functions [U [22j , where the GLRMQ terms are 
present. This allows us to test the dynamic responsible 
for the rise of total pp{p) cross sections with energy in 
the presence of saturation effects. 

In order to compare these two different regimes, we use 
here GRV98(LO) [1| and CTEQ6L Q for the parton 
densities, as a baseline. They are governed by DGLAP 
evolution equations (linear regime), where, therefore, sat- 
uration effects are clearly absent. 

At this point, we would like to call attention for the lat- 
est LHC results: they have provided very valuable infor- 
mation on high-energy multiparticle production, improv- 
ing our theoretical understanding of soft and semi-hard 
parton dynamics, and showed the need for adjustments 
and even the reformulation of hypotheses employed in 
models that propose to establish the behavior of hadronic 
cross sections with energy. 

It is also important to note that the recent cosmic ray 
data from HiRes and P. Auger collaborations have al- 
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lowed a deeper understanding about the nature of pro- 
duced particles at very high energies and stimulated 
many discussions between the accelerator and cosmic-ray 
communities on common issues in these areas [IH, [25j . 

The energy dependence of total hadronic cross-sections 
is probably the most important question for the cosmic- 
ray community. The phase space regions of relevance to 
the development of air showers are not directly accessible 
in the currently accelerator experiments, and, because of 
that, descriptions and interpretations of data in cosmic 
rays physics at high energies depend crucially on the pre- 
dictions coming from phenomenological models. 

For this reason it is very interesting to test the range of 
the model presented in this work and to verify if it per- 
mits a satisfactory description of data to other processes. 
Obviously, our first choice is the inelastic proton-air cross 
section. 

In the Glauber multiple collision model [26| the inelas- 
tic proton- nucleus cross-section, cr^, ; (s), can be derived 
in the eikonal limit (straight line trajectories of collid- 
ing nucleons) from the corresponding inelastic nucleon- 
nucleon (NN) cross-section, ufj^^s). At the center-of- 
mass energy t/s and the geometry of the pA collision, 
it is simply determined by the impact parameter of the 
reaction: 



pA 
inel 



d 2 b 



1 



(5) 



The thickness function 7^(6) (= J dzpA(b, z), with 
J d 2 bTA(b) = A)gives the number of nucleons in the 
nucleus A per unit area along a direction z, separated 
from the center of the nucleus by an impact parame- 
ter b, and paQ>,z) is the nuclear density of nucleus A. 
In what follows we have used A a i r = 14.5, pA(b, z) = 
po {1 + exp [(r - i? A )/ao]} _1 with r = Vb 2 + z 2 , R A 
1.1?) A 1 /* - 1.61 A- 1 / 3 (fm), a = 0.75 fm. 

Our second choice to test the model are the jp and 77 
cross sections. The data on these processes can be de- 
rived from the pp forward scattering amplitude using Vec- 
tor Meson Dominance (VMD) and the Additive Quark 

Model with the introduction of a probability {P^ad 1 ^) 
that the photon interacts as a hadron 0, [2(| ■ 

Assuming, in the spirit of the VMD, that at high en- 
ergies the photon behaves as a hadronic state composed 
by two quarks, after the substitutions a 



s,h 



%a s ' h and 



Ms, ft "~ ^ \J 5A*s,h m both, soft and hard components of eq. 
©, the jp cross section can be written as: 



aZ{s)=2PZ d j d'b{l-e 
x cos [Rex lp (b,s)}} , 



-Imx 1P {b,s) 



(6) 



4 



1 60 r {■ I 1 — I ' — I 1 — I 1 — I / M [ ■ I ' — I 1 — I ' — I I ' — n™ 




FIG. 1. Total pp and pp cross sections. Solid and dashed lines represent the results of the model (eq. JT])) and the experimental 
data are from references [27l-[3l|]. 



To obtain the 77 cross section we apply the same pro- 
cedure as before making the substitutions 

a a,h ± a s,h 

and (i a>h ->• 

= J d 2 b{l e -^x-(M) 

x cos[Re X ^(b, S )}}, (7) 

where = (J™ )>. 

The simplest VMD formulation only assumes the light- 
est vector mesons and in this case P^d * s gi ven by 

P had = T,V=p,u,$ ^jf^' Where "e™ iS the Q ED Coupling 

constant and f v is the ^y — V coupling. In this work, how- 
ever, we consider P^ d as a free parameter fixed by the 
low energy 7p data. 

III. RESULTS AND DISCUSSIONS 

First of all we would like to emphasize that the PDFs 
used here have different kinematic limits. Their ranges 
of validity are: 1CT 5 < x < 1 and 1.4 GeV 2 < Q 2 < 
10 4 GeV 2 (EHKQS); 1CT 9 < x < 1 and 0.8 GeV 2 < Q 2 < 
10 6 GeV 2 (GRV98); 10~ 6 < x < 1 and 1.69 GeV 2 < 
Q 2 < 10 s GeV 2 (CTEQ6L). Although the last distribu- 
tion incorporates some physics of small- x , only the first 
one is governed by nonlinear effects of QCD. Second, in 
our calculations were considered only contributions which 
respect the kinematic ranges (for both, x and Q 2 ) of the 
distributions above. 

The strategy adopted in our calculations is as follows: 
while the soft contribution of the model is parameter- 
ized to describe the pp and pp scattering at low energies, 
at higher energies the hard parameters are fixed in the 
nonlinear regime (using the EHKQS parton distribu- 
tion functions) to fit the latest LHC data and the re- 



sults from HiRes and P. Auger collaborations (converted 
data) . 

In all figures bellow, we have adopted for left pan- 
els the values p\ min — lAGeV 2 , [i gg = 1.65 GeV and 
[i qq = 1.6 GeV and, for right panels, the values p\ . = 
1.69 GeF 2 , fj, gg = 1.6 GeV and n qq = 1.5 GeV . In" all 
cases we have fixed ^ 2 so f t — 0.71 GeV 2 (see eq. ([3])) for 
the soft component. Solid curves show the results ob- 
tained with nonlinear evolution and dashed curves rep- 
resent the results obtained with these same parameters 
(fixed considering saturation effects) in the linear regime. 
This procedure allows us to compare directly the two dif- 
ferent dynamics used in our studies. 

The relevant experimental data showed in the next fig- 
ures considers only the most quoted ones in the literature 
and can be found in the references [2(| [27l - l35l | . 

Figure Q] shows the results of pp and pp total cross sec- 
tions from the present model (eq. ([T])). As can be seen 
the results obtained with nonlinear evolution produce 
a good description of data [27H3l| . It should be noted 
that, with linear evolution, a better description of data 
could be obtained with different choices of parameters 
above. However, in this case, the comparison between 
the nonlinear and linear regimes, would make no sense. 
Comparing the solid curves with the dashed ones, we 
conclude that the rise of these cross sections at high en- 
ergies is still dictated by the growth of gluon densities in 
small- a; but it is suppressed due to the recombination of 
gluons, leading to a softening of these cross sections. 

Figure [2] shows the ratios between the results pre- 
sented in Figure [1] and the parametrization er(s) = 
30 + 0.22 ln 2 (s/s ) (mb), with s = lGeV 2 , similar 
to those shown in Particle Data Group (best fit) [27| . 
We have adopted it here because its "universal term", 
cx ln 2 (s/ So), permits a more satisfactory description of 
data at high energies, including the latest LHC data and 
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FIG. 2. Ratios between total pp and pp cross sections (solid and dashed lines of Figure [T]) and the parametrization presented 
in the text. 
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FIG. 3. Total inelastic pp and pp cross sections. Solid and dashed lines represent the results of the model (eq. (0) and the 
experimental data are from references [27j - [3T| ] . 



the results from HiRes and P. Auger collaborations (con- 
verted data). This simple exercise, however, indicates us 
a very interesting result: the nonlinear dynamics favor 
more effectively than the linear one the respect to the 
Froissart's limit by the hadronic cross sections. 

Our results for the inelastic pp and pp cross sections 
(eq. ([T])) are showed in Figure [3] and compared to the 
experimental data, which include the most recent re- 
sults from P. Auger Collaboration [3l[ and LHC fpj . 
and the oldest ones from the Particle Data Group [27| 

ip&{") = "Sf^OO - ^\s)). Evidently, the de- 
scription of data is not so good as before. The relatively 
weak growth observed in these cross sections is related, 
on the one hand, to the kinematic limits of PDFs (which 
exclude lower values of PT mi „ and x than those used here) 



and, on the other hand, because, as in other models (ref- 
erences cited in this work), the imaginary part of our 
eikonal function, which controls these cross section, does 
not contain diffractive processes. However, albeit com- 
plemental, this figure shows the role played by the satu- 
ration effects in describing these experimental data, and 
tells us that the inclusion of diffractive processes is really 
needed. 

Figure @] shows our results for the inelastic proton-air 
cross sections given by eq. ([5]). They take into account 
only the geometry of the collision p — Air and the semi- 
hard dynamics of the model, via the imaginary part of the 
eikonal function, contained in the inelastic proton-proton 
cross sections. As discussed above, the parameters of the 
model were fixed only through the total experimental pp 
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FIG. 4. Inelastic proton-air cross sections. Solid and dashed lines represent the results of the model (eq. ([5])) and the 
experimental data are from references [2p|. l3ll. [3^] . 
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FIG. 5. Total gamma-proton cross sections. Solid and dashed lines represent the results of the model (eq. ([6])) and the 
experimental data are from references [27l.l34l]. 



and pp data cross sections, which do not provide a so 
satisfactory description of inelastic pp collisions, because 
diffractive processes are absent and we have restricted 
ourselves to the kinematic limits of the EHKQS distri- 
bution functions in the nonlinear regime. Despite that, 
the results obtained with nonlinear evolution describe 
satisfactorily the experimental data [2(J HH, |33[ , in par- 
ticular from HiRes and P. Auger collaborations. We con- 
sider their as a lower limit of our calculations, but, in our 
opinion, they may be compared, quantitatively, to the 
predictions of other models (see discussions presented in 

Figures [5] and H] show our results for the total 7p and 



77 cross sections, respectively given by equations ([5]) and 
([7]). The probability that in these collisions the pho- 
ton interacts as a hadron was fixed in (P^d ~ 1/221). 
The experimental data are from references [271 [34| and 
[35j , respectively. The good agreement between our solid 
curves (nonlinear regime of the model) and the data is 
a nice surprise considering that these cross sections are 
determined by the same parameters used in pp and pp 
collisions and by the changes coming from the "weights" 
introduced by the VMD and Additive Quark Model in 
the present version of the eikonalized minijet model. 
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IV. CONCLUSIONS 

In this work we have used an eikonalized minijet model, 
where the effects of the first nonlinear corrections to the 
DGLAP equations are taken into account, to investigate 
the behavior of total hadronic and photonic cross sec- 
tions. 

Firstly, it is important to note that the same approach 
was used to investigate the behavior of cross sections 
measured in accelerators and cosmic rays experiments. 
Secondly, we would like to call attention for two very 
interesting and important results of our studies: a) the 
introduction of nonlinear corrections to the model allow 
for a very satisfactory description of total hadronic and 
photonic cross sections, simultaneously; b) even though 



small and only manifest themselves at extremely high en- 
ergies, the saturation effects attenuate the strong growth 
of total hadronic cross sections obtained by conventional 
eikonal models and contribute significantly to an asymp- 
totic behavior which satisfies more comfortably the Frois- 
sart bound. 

Finally, the fact that we have achieved a reasonable 
success in using this model, encourages us to test it (work 
in progress) as the description of other observables. 
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